Orbital order in bilayer graphene at filling factor v = —1 



O 
(N 



(N 



O 

o 



> 

ON 

\o 

rn 
O 

o 



X 



R. Cote/ Jules Lambert,^ Yafis Barlas,^ and A. H. MacDonald^ 

^ Departement de physique, Universite de Sherbrooke, Sherbrooke, Quebec, J IK 2R1, Canada 
^National High Magnetic Field Laboratory and Department of Physics, 
The Florida State University, Florida 32306, USA 
''Department of Physics, The University of Texas at Austin, Austin, Texas 78712, USA 

(Dated: March 31, 2010) 

In a graphene bilayer with Bernal stacking both n = Q and n — 1 orbital Landau levels have 
zero kinetic energy. An electronic state in the A*' = Landau level consequently has three quantum 
numbers in addition to its guiding center label: its spin, its valley index K ox K' , and an orbital 
quantum number n = 0, 1. The two-dimensional electron gas (2DEG) in the bilayer supports a wide 
variety of broken-symmetry states in which the pseudospins associated these three quantum numbers 
order in a manner that is dependent on both filling factor v and the electric potential difference 
between the layers. In this paper, we study the case oi v = —1 m an external field strong enough to 
freeze electronic spins. We show that an electric potential difference between layers drives a series of 
transitions, starting from interlayer-coherent states (ICS) at small potentials and leading to orbitally 
coherent states (OCS) that are polarized in a single layer. Orbital pseudospins carry electric dipoles 
with orientations that are ordered in the OCS and have Dzyaloshinskii-Moriya interactions that can 
lead to spiral instabilities. We show that the microwave absorption spectra of ICSs, OCSs, and the 
mixed states that occur at intermediate potentials are sharply distinct. 

PACS numbers: 73.21.-b,73.22.Gk,78.70.Gq 



I. INTRODUCTION 

Semiconductor double-quantum-well systems in a 
quantizing magnetic field develop spontaneous inter- 
layer coherence when the wells are brought into close 
proximityji Spontaneous coherence leads to a variety of 
fascinating transport effects including counterflow super- 
fluidity and anomalous interlayer tunneling, and to un- 
usual charged excitations such as merons. A convenient 
way to describe these ground states is to use a pseudospin 
language in which the which layer degree-of-freedom is 
mapped to a 5 = 1/2 pseudospin. In this language, the 
ground state of a bilayer at total filling factor v — 1 is 
an easy-plane pseudospin ferromagnet. At higher filling 
factors, still more exotic states occur , for example states 
in which the pseudospin orientation varies in space and 
a charge-density-wave is formed.^ 

Interest has recently been growing in the strong- 
magnetic-field ordered states of graphene bilayers. Single 
layer graphene^ is a two-dimensional honeycomb lattice 
network of carbon atoms. Bilayer graphene^""— consists of 
two graphene layers separated by a fraction of a nanome- 
ter. In the normal Bernal stacking structure, one of the 
two honeycomb sublattice sites in each layer has a near- 
neighbor in the other layer, and one does not. This ar- 
rangement produces^ a set of Landau levels with ener- 
gies En = ±fihjl^\N\ {\N\ + 1) where uj*^ is the effective 
cyclotron frequency and N = 0,±1,±2,.... All Landau 
levels except = are four-fold degenerate; electronic 
states are specified by N , valley-index [K or K') and 
spin-index, in addition to the usual label used to specify 
guiding center states within a Landau level. 

The A^ = Landau level has an additional two- 
valued quantum degree-of-freedom because states with 



both n = and n ~ 1 Landau-level character have zero 
kinetic energy. Most of the new physics discussed in this 
paper is related to the property*" that electric dipoles 
can be constructed by forming wavefunctions with co- 
herence between n — Q and n — \ components. A sec- 
ond peculiarity of the A^ = state is that wavefunctions 
associated with the K valley are localized in one layer, 
while wavefunctions associated with the K' valley are lo- 
calized in the opposite layer. Layer and valley indices 
are thus equivalent. It is convenient to use pseudospins 
to represent both layer (or equivalently valley) and the 
Landau-level orbital character degrees of freedom. The 
wavefunction for an electron in the A^ = Landau level is 
therefore the direct product of a standard guiding center 
factor and three spinors that capture its dependence on 
spin, layer, and orbital-Landau-index (n) character. We 
refer to the final spinor as the orbital spinor, and to the 
set of eight Landau levels with zero kinetic energy as the 
bilayer graphene octet In neutral graphene the octet is 
half-filled at all magnetic field strengths. 

The presence of the octet in bilayers is revealed exper- 
imentally by a jump in the quantized Hall conductivity^ 
from — 4(e^//i) to 4(e^/ft.) when the charge density is 
tuned across neutrality in moderately disordered sam- 
ples. In a recent paper— some of us predicted that quan- 
tum Hall effects would appear at all integer filling fac- 
tors between v = —4 and v — A m samples of quality 
sufficient^ to make interactions dominant relative to un- 
intended disorder. Electron-electron interactions acting 
alone are expected to lift the degeneracy of the bilayer 
octet and induce gaps at the Fermi level by producing a 
set of spontaneously broken symmetry states with spin, 
valley and orbital pseudospin polarizations. The octet 
degeneracy lifting is expected^ to follow a set of Hund's 
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rules in which spin polarization is maximized first, then 
layer polarization to the greatest extent possible, and fi- 
nally orbital polarization to the extent allowed by the 
first two rules. Hall plateaus at all integer filling factors 
intermediate between v = —4: and v — -i have indeed now 
been discovered in experimental studies of suspended bi- 
layer graphene samples and bilayer graphene on Si02/Si 
substrates,-^? opening up the opportunity to study a rich 
and still relatively unexplore d^ ^'^^ family of novel bro- 
ken symmetry states. The odd filling factor cases are 
expected to be most interesting because all three pseu- 
dospins are expected to be polarized. The present paper 
focuses on the physics associated with the competition 
between layer and orbital pseudospins at Landau levels 
1/ = — 1 and J/ = 3 at field strengths sufficient to produce 
maximal spin polarization and reduce the importance of 
Landau-level mixing. In this limit a negative filling factor 
h' is equivalent to a positive filling factor v -\- A since the 
two states differ only through the presence in the latter 
case of inert filled majority spin Landau levels. 

In a previous paper^A we studied the quantum Hall 
states which occur at j/ = —3 and v = 1 'm the same 
strong field regime, emphasizing the key role played by 
the potential energy difference between graphene layers 
which we refer to here as the bias potential A^. The 
V = —i ground state at zero bias is an inter- layer coher- 
ent state with orbital index n = that supports counter- 
flow superfluidity. One particularly interesting property 
of this state is that the superfluid density, the coefficient 
that relates the counterfiow supercurrent to the spatial 
gradient of interlayer phase, vanishes. Correspondingly, 
the state's Goldstone mode dispersion is quadratic in 
wa vector q, in contrast to the linear dispersion found in 
coherent semiconductor bilayers and in standard super- 
fiuids. We also found that the uniform ground state has 
a long-wavelength instability at any non-zero potential- 
difference bias As < A^'' where A^-* is the critical bias 
at which all iV = charge is tranferred to a single layer. 
In Ref. [H we argued that the instability is probably 
towards a state in which the direction of the inter-layer 
pseudospin varies in space. For larger bias As > A^'^ the 
ground state is uninteresting; the charge is completely 
in one layer (or valley) and in the orbital state n — Q. 
The orbital pseudospinwave mode corresponding to tran- 
sitions between the n = and n ~ 1 orbital states is 
gapped at a frequency w = fia;*As/7i where 71 is the 
inter-layer tunneling energy in the Bernal stacking. This 
mode, which is an intra-Landau level excitation, has a 
finite oscillator strength and will absorb^^ electromag- 
netic radiation. This behavior contrasts with the stan- 
dard Kohn's theorentii behavior in normal 2DEG's which 
implies that only inter-Landau level excitations produce 
absorption. 

Surprisingly the phase diagrams for v = ~1 and v ^ 3 
states differ qualitatively from the corresponding 1/ = —3 
and z/ = 1 phase diagrams. The source of the difference 
is a competition in the v = —1 case between interac- 
tion and single-particle effects which are reinforcing in 



the ly = —3 case. The end result is that the large A^ 
ground state at v = —1 places electrons in a coherent 
combination of n = and n — 1 orbital states, and that 
electric dipoles are consequently spontaneously present 
in the ground state. This paper analyzes the dependence 
of bilayer properties on As and explores some of the con- 
sequences of the unusual orbitally ordered dipole state. 

At small bias we find that the bilayer 's 1/ — —1 ground 
state is an inter-layer coherent state, much like the cor- 
responding v = —3 state except that the coherence is 
between orbitals with n = I character. This state has 
a gapless pseudospin wave mode with linear dispersion, 
like coherent semiconductor bilayers. The state also has a 
gapped orbital pseudospin collective mode. Because the 
orbital spinor carries an electric dipole, this mode has 
a finite oscillator strength and absorbs electromagnetic 
radiation, again much like the v = —3 case. This mode 
should be visible in a microwave spectoscopy experiment. 

Inter-layer coherence decreases with bias until a new 
ground state is reached that has both inter-layer and or- 
bital coherences. In this mixed state, the low-energy or- 
bital and inter-layer pseudospin modes are both gapped. 
Because the modes are coupled, both show up in the mi- 
crowave absorption spectra. The collective excitations 
are highly anisotropic in this phase and we find that the 
intensity of the absorption depends strongly on the ori- 
entation of the electric field of the incident microwaves. 

The new physics of the 1/ ~ —I case emerges in its 
simplest form at still stronger bias potentials. Both or- 
bital levels in the bottom layer are then completely filled 
while only one of the two top layer Landau levels is filled. 
Spontaneous orbital coherence then develops in the top 
layer. This spontaneous orbital coherence leads to a gap- 
less orbital pseudospin mode. Some of the properties of 
this state have been studied independently in a recent 
paper by Shizuya-^, who also pointed out that orbital 
coherence is responsible for the existence of a finite den- 
sity of electrical dipoles with a net polarization. These 
dipoles collectively and spontaneously point in some ar- 
bitrary direction in the x — y plane. As discussed by 
Shizuyaii^ their orientation can however be controlled by 
an external electric field parallel to the plane of the bi- 
layer. In this paper, we use an effective pseudospin model 
to highlight other interesting features of the orbitally- 
coherent state. In particular we demonstrate the pres- 
ence of a Dzyaloshinskii-Moriya (DM) interaction-'^^ be- 
tween orbital pseudospins and show that it leads to an 
anisotropic softening of the orbital pseudospin mode at 
a finite wavevector. For strong enough inter-layer bias, 
the DM induces an instability toward a pseudospin spi- 
ral state. The orbital pseudospin mode in the high bias 
regime is gapless and will lead, in the presence of dis- 
order, to strong absorption of electromagnetic waves at 
very small frequencies. 

Our paper is organized in the following way. In Section 
II, we discuss the non-interacting states of the graphene 
bilayer within a two-band low-energy model. Here we in- 
troduce the aspect of the electronic structure that is re- 
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sponsible for interaction and band effects which are com- 
peting aX v = —\ and are reinforcing aX v = —3. In Sec- 
tion III, we derive the Hamihonian of the graphene two- 
dimensional electron gas (2DEG) truncated to = 
levels in the Hartree-Fock approximation. We use this 
Hamiltonian to derive the equation of motion for the 
single-particle Green's function in Section IV and to ob- 
tain the order parameters for the various phases which 
occur at v — —1. Section V describes the generalized 
random-phase approximation (GRPA) (or equivalently 
the time-dependent Hartree-Fock approximation (TD- 
HFA)) that we use to derive the collective excitations. 
The phase diagram for = — 1 as a function of bias is 
obtained in Section VI. We then study the collective ex- 
citations of the inter-layer coherent phase in Section VII 
and those of the orbital coherent phase in Section VIII. 
Finally microwave absorption in the different phases is 
studied in Section IX and we conclude with a brief sum- 
mary and some suggestions for future work in Section 
X. 



II. EFFECTIVE TWO-BAND HAMILTONIAN 



In a magnetic field, 



Him 



-e^As + e/3As (1 + a^a) 



(2) 

where we have defined the orbital ladder operators a — 
e {P^ - iPy) /V2h, at = f (F^ + iPy) with the mag- 

netic length I — ^Jhc/eB and the parameter j3 — 
nw*J"fi = 6.144 X 10-3B(Tesla). The effective cyclotron 
frequency w* — eBjm'c. At zero bias, the Landau 
levels have energies E% = ±ihujl^\N\ {\N\ + 1) with 
A^ = 0,±1,±2,... 

In this paper, we study the phase diagram of the 2DEG 
in the A^ = Landau level. While levels with |A^| > are 
four-fold degenerate (counting spin and valley quantum 
numbers), level A'^ = has an extra orbital degeneracy 
due to the fact that Landau level orbitals n — and n = 1 
have zero kinetic energy. The states in A^ = are thus 
member of an octet of Landau levels that are degener- 
ate if we neglect the Zeeman and bias potential energies. 
We assume that the Zeeman coupling is strong enough 
to assure maximal spin-polarization, which allows this 
degree-of-freedom to be neglected. The eigenfunctions 
and corresponding energies for N — are then given by 



In a graphene bilayer with Bernal stacking, the two 
basis atoms of the top layer are denoted by Ai and Bi 
and those of the bottom layer by A2 and B2 with atoms 
Ai sitting directly above atoms -62- The band structure 
of the bilayer is calculated using a tight-binding model 
with in-plane nearest-neighbor tunneling (with strength 
7o — 2.5 eV) and Ai — B2 tunneling (with strength 71 = 
0.4 eV). (For a review of bilayer graphene, see Ref. [l6|.) 
The low-energy {E << 71) excitations of this model for 
electrons in the valleys with K — (— 1,0) and 
^ = ±1 can be studied by using the effective two-band 
model developped in Ref. ,5;. Using the basis {|^2) , l-Bi)} 
for and {|_Bi) , IA2)} for H°j^, the effective two-band 
Hamiltonian derived in Ref. H] is 



H^K - 



_f Ab , £Ab^ 



71 2m* 



P+ 




where p± — Px^ipy and p = —ilW. In this equation, A^ 
is the bias potential between the two layers, the effective 
mass m* = 2fi^7i/37oaQ = 0.054mo with mo the bare 
electronic mass, ao = Vic is the triangular lattice con- 
stant and c — 1.42 A is the distance between neighboring 
carbon atoms in the same plane. The kets IA2) , \Bi) 
correspond to the atomic sites in different layers that are 
not directly above one another. 

The Hamiltonian of the 2DEG in a perpendicular 
magnetic field is obtained by making the substitution 
p ^ p -t- eA/c = P/fi (with e > 0) in Eq. (P. The vec- 
tor potential A is defined such that V x A = B =Bz. 







ho.x (r) 


hi,x (r) ; ' 
for the K valley and by 



(3) 



^Ab-Z^Ab, (4) 



ho.x (r) 


hi,x (r) 




K',n=0,X 



K',n=l.X 



Ab, (5) 
Ab+/3Ab (6) 



for the K' valley, using this time the basis {IA2) , I -Si)} 
for all states. It is quite clear from these equations 
that the valley K(K') eigenstates are localized in the 
top(bottom) layer. For A^ = 0, the layer index is thus 
equivalent to the valley index. (For |A^| > 0, the spinors 
have different orbital indices n in different layers.) The 
functions hn,x (r) = er^^^^^ ipn [x — X) /y^X^ are the 
Landau gauge (A = (0, Bx, 0)) eigenstates of an electron 
with guiding center X, and (/?„ (x) is the wave function of 
a one-dimensionnal harmonic oscillator. Note that with 
our choice of gauge, the action of the ladder operators on 
the states tfn {x) is given by a^(fn {x) = l^Pn+i {x) 

and aipn (x) = -iy/n(pn-i {x) . 

At finite bias, the parameter /3 << 1 lifts the degener- 
acy between the two orbital states as we show in Fig. 
[T] The splitting is however very small. For positive 
bias, the n — orbital state in the bottom(top) layer is 
lower(higher) in energy than the n = 1 orbital state. The 
orbital states n = 0, 1 form a two- level system in each 
valley and we associate them with an orbital pseudospin. 
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Similarly, the two states ±K are associated with a val- 
ley pseudospin. We remark that the effective two-band 
model slightly overestimates the gap A = — i?n=o 

between the n — and n — 1 orbital states. The inset in 
Fig. [T]shows the difference between the gap calculated in 
the two-band model and in the original four-band system. 
In the region where the DM interaction driven instabil- 
ity occurs, the difference between the two gaps is however 
very small. 



and 



*K' (r) - 



X 



E 

X 



ho,x (r) 


It-i,x (r) 




CK'.Xfl 



CK'X.I, 



(8) 



so that the Hartree-Fock Hamiltonian is given by (here 
and in the rest of this paper, we use the convention that 
repeated indices are summed over) 
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FIG. 1: Non-interacting energy levels in Landau level A'^ = 
0. The inset shows the behavior of the gap evaluated in the 
original (solid line) and two-band (dashed line) models with 
bias. 



HhF — N<pEa,nPn% (0) 



(9) 



(q) {Pn\\, (-q)) p'nlu (q) 



^/ ^ 711,714,713,712 

q 

where is the Landau level degeneracy and all energies 
are measured in units of j where k is the effective di- 
electric constant at the position of the graphene layers. 
The single-particle energies Ea,n include capacitive con- 
tributions and are defined by 



V d 
2I 



(10) 



with Va the number of filled levels in valley v = v + A 
the total number of filled levels, a, & = ±1 the valley (or 
equivalently layer) index and n = 0, 1 for the two orbital 
state indices. In deriving Eq. ([S]), we have taken into 
account a neutralizing positive background so that the 
q = contribution is absent in the Hartree term. This 
convention is indicated by the bar over the summation. 
Note that for positive bias, the bottom layer (K' valley) 
is at a lower potential than the top layer (K valley). 
The density operators in Eq. are defined by 



III. HARTREE-FOCK HAMILTONIAN 



We now add the Coulomb interaction to the non- 
interacting Hamiltonian H%^. We assume that the mag- 
netic field is strong enough so that we can restrict the 
Hilbert space to the A*" = Landau level and neglect 
Landau level mixing. We also assume the 2DEG to be 
fully spin polarized (we comment on this later). We write 
the electron field operator as 



- E ( (r) ) 



(7) 



(r) 



ck,x,i 



<,77. (q) = ]^ E 



^-iq^{Xi+X2) 



(11) 



^1 ,^2 



^cl,Xi,7iiCb,X2,7i2'5xi,X2+g^f2, 



where c^^ Xx ni creates an electron in state {a, Xi,ni) in 
the Landau gauge. The intralayer {H,X = H"^" , X 
and inter-layer (^H,X ^ H'^^^.X"^^^ Hartree and Fock 
interactions are given by 
1 



7ii, 712, 713, 774 (q) 



Km,n2 (q) -^713,714 (-q) , (12) 



^711,712,713,714 (q) 

and 

Hni,n2,n3,n4 (q) 

-X^7ii,n2,7i3,7i4 (q) 



J ^771,712,713,714 (p) ^ ^ (l^) 



-f^7li,7l2,7l3,7l4 (q) 6 
1^2 



— qd 



(14) 



J ^"1. "2, 713, 714 (P) S ^ il^) 
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where d = 3.337 A is the inter-layer separation in the 
Bcrnal stacking. The form factors which appear here, 



i^o.o (q) 
Kui (q) 
Ki,o (q) 
^04 (q) 



exp 
exp 



-q^e 

4 



q^e 



(qy + iqx)i \ I -q 

[-qy + iqx)(. 
V2 



exp 



4 

-qH^ 



(16) 
(17) 
(18) 
(19) 



capture the character of the two different orbital states. 
Detailed expressions for the Hartree and Fock interac- 
tions parameters are given in Appendix A. 



IV. ORDER PARAMETERS AT INTEGER 
FILLINGS 

The states with no pseudospin texture at inte- 
ger filling factors have uniform electronic density and 
density-matrices that vanish for q 7^ 0. Letting 
(Pnf.n. (q = 0)) ^ (p;^f,„J , the Hamiltonian of Eq. © 
reduces to 



H 



HF 



a. a 



(20) 



(0) 



These order parameters are conveniently calculated by 
defining the time-ordered Matsubara Green's function 



G^n. {r) = -^Y. (TrCa.,n.x (t) 4.„„x (0)) , (21) 



since, at time zero, we have 



(22) 



In the Hartree-Fock approximation, the equation of 
motion for the single-particle Green's function is 



IHiUJr, 



{Ea^n, - ^^)] G'^tn. K) (23) 
where w„ is a fermionic Matsubara frequency and 



K:k.,n,n'io){f^a)- (24) 



jja.c 



The system of Eqs. (|23l) can be solved in an itera- 
tive way by using some initial values for the parameters 
{{Pn2,ni)} ■ I^^f- B solved this equation keeping 
valley, orbital, and spin indices. We showed that the so- 
lutions of the Hartree-Fock equations for the balanced 
bilayer (A^ = 0) follow a Hund's rules behavior. The 



spin polarization is maximized first, then the layer po- 
larization is maximized to the greatest extent possible, 
and finally the orbital polarization is maximized to the 
extent allowed by the first two rules. In the absence of 
bias, the ordering of the first four states (with spin up) 
is given by 



1^,0) 
1 5,1) 
\AS, 0) 
\AS, 1) 



^|i^,o) + -L|^',o) 



V2 



\K,l) 



V2 



\K'A) 



(25) 
(26) 
(27) 
(28) 



in this order. The next four states follow the same order 
but with spin down. The occupation of these eight states 
are given by the filling factor v ranging from v = —3 
(state I S, 0) with spin up fully filled) to = -1-4 (all eight 
states filled). 

To simplify the notation, we define 

\K,0) ^ |1), (29) 

\K,l) ^ |2), 
\K',0) ^ |3), 

I^M) ^ |4), 

so that ^ {p^,J) with ij = 1,2,3,4. 

It is easy, using Eq. ([53]), to prove the sum rules 



where 



'^i = {Pi,i) 



(30) 



(31) 



V. COLLECTIVE MODES IN THE 
GENERALIZED RANDOM-PHASE 
APPROXIMATION 

In order to compute the collective excitations, we de- 
fine the two-particle Matsubara Green's function 



AtIi ,Tl2 ,714 / 



(32) 



where again rii =0,1 are orbital indices and a,b,c,d 
are valley indices. To derive the equation of motion 
for these response functions in the Generalized Random- 
Phase Approximation (GRPA), we proceed in the fol- 
lowing way. We first derive the equation of motion for 
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X in the Hartree-Fock Approximation (HFA) using the 
Heisenberg equation of motion 



[H - mTV, (. . .)] . 



(33) 



where H is the Hamiltonian of Eq. ([9]) with the av- 
erages removed, /i is the chemical potential and N the 
number operator (not to be confused with the Lan- 
dau level index). After evaluating the commutators, 
we linearize the resulting equation by writing p (q) — > 
{p {<\))fjpA + (q)- We get the GRPA equations of mo- 
tion by keeping the terms up to linear order in 5p (q) . In 
the homogeneous states at integer filhngs, {p{q))^p = 
(p il)) HF ^<i,o, so that we get the set of equations 

= ^ {PntnJ 5n^,njb,c - ^ (p^3 ) '^a,d'5„i ,„4 (34) 

, , (0) / p"^!'' , ) x^T' (q, ao 



(0) A 



a .b 



(0)a,a',c,d , ^ . 
Aril, 71^, 713,714 V^i "n^ ) 



and 



(35) 



a,b,c,d / .Q \ 

A711 ,712,713,714 \Hi ^'n^ 

_ iO)a,b,c,d / .Q N 

A7ll,7l2,7l3,7t4 \4l '^'71^ 

i}.-yiO)a,b,e,e / . q \ fje.g /^ 

^ ^A7li,n2,7l5,7l6 VMi ^'71^ -'-'715,715,717,713 \HJ 

g,g,c,d / Q ^ 
^ An7,7is,7i3,n4 V4)^'n; 

^A7ii,7i2,7i5,7i6 ^ ) ^ Us ,na ,717 ,ne VM^ 

^Xn'7,7lS,7l3,7l4 (qi ^71 ) , 

where f7„ is a bosonic Matsubara frequency. The re- 
tarded response functions are obtained, as usual, by tak- 
ing the analytic continuation if2„ cj + iS. 

By defining super indices A, B — 1,2, 3, ...,16 repre- 
senting the combinations (a, ni; 6, 77.2), (c, ng; d, 714) , etc., 
we can represent the response functions and interactions 
matrices as 16 x 16 matrices and then write the GRPA 
equation in the matrix form: 



identified in this way, three are Pauli-blocked at v = —1 
and appear in our calculations as dispersionless modes 
that have zero weight in all physical response properties. 
The three remaining excitation modes are physical and 
for V = —1 correspond to the interaction-coupled transi- 
tions indicated in Fig. [51 Note that, in the limit q — > 00, 
H'''3,X'='f so that Eq. gives x ^ X°- In this 

limit, the collective mode frequencies correspond to tran- 
sitions between eigenstates of the Hartree-Fock Hamilto- 
nian Hhf as illustrated in Fig. [5] for the case v — —1 
and zero bias. The energy of these eigenstates include the 
non-interacting energies and the self-energy corrections. 
The w (q = 0) limit of these sames modes, however, also 
includes the polarization and excitonic corrections that, 
in a Feynman diagram description of the GRPA, are cap- 
tured by bubble and ladder diagram summations. These 
effects make the modes dispersive. 
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FIG. 2: At f = —1, the first three energy states of the 
Hartree-Fock Hamiltonian are filled and the only possible 
transitions are those indicated by the three arrows. The fre- 
quency of these transitions corresponds to the frequency of 
the collective excitations calculated in the GRPA in the limit 
in which the wavevector g — >■ cxd. We illustrate here the case 
of zero bias. 



[(u; + z(5)/-F(q)]x(q,c^)-S(q), (36) 

where B,I,F,x are 16 x 16 matrices (with / the unit 
matrix) . The matrices F (q) and B (q) depend on the 

(Pn° 711 evaluated in the HFA. We will give later the 
precise form of these matrices for the phases studied in 
this paper. 

The frequencies of the collective excitations are given 
by the eigenvalues of the matrix F (q) . There are in 
total 4 zero modes, corresponding to unphysical intra- 
level transitions, and 12 non-zero modes corresponding to 
inter-level transitions. The latter occur in six positive- 
negative energy pairs, corresponding to excitation and 
deexcitation partners. Of the six collective excitations 



VI. PHASE DIAGRAM AT FILLING FACTOR 

The properties of the ground state at = —3 have been 
studied in detail in Ref. [U At zero bias, the ground 
state has all electrons in the 15, 0) state and can be de- 
scribed as an XY layer-pseudospin ferromagnet with or- 
bital character n = 0. The bias acts as an effective 
external magnetic field that forces the layer pseudospins 

(c) 

out of the x — y plane. Above a critical bias , all elec- 
trons are in the bottom layer and the layer-pseudospin 
is correspondingly fully polarized. The ground state is 
then given by \K' , 0) and is unchanged if the bias is fur- 
ther increased. The v = —1 state on which we focus 
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here differs from tfie v = —3 state because of a com- 
petition between single-particle and interaction energy 
effects which emerges only in the former case. As il- 
lustrated in Fig. (HI, single-particle effects captured by 
the bilayer effective Hamiltonian favor occupation of the 
n — 1 orbital when three of the octet's eight levels are 
occupied (v = —1). Note that this tendency is indepen- 
dent of the sign of A^. Exchange interactions, on the 
other hand, always favor a state in which as many n = 
orbitals as possible are occupied. As we explain below, a 
compromise is reached by forming a state with coherence 
between n — and n — 1 orbitals. This physics is en- 
riched by the same tendency toward interlayer coherence 
which occurs at v = —3 and^ in semiconductor bilayers. 
Indeed, our calculations show that the phase diagram at 
v = —1 is much more complex than at j/ = — 3. 



A. Inter-layer-coherent state 

At zero bias, numerical solution of the HFA equations 
leads to occupied \S,0) , \S, 1) , and \AS,0) states. The 
order parameters are then given by 



The energy of this inter-layer-coherent state (ICS) for 



As < A^j^' is 



^ = ^A^ (1 - 2/3) (P,i) 



(42) 



-lx,{0)-^x^ (0)-^xi6 (0) 
+ ^(^-^16 (0)) 

-^xie (0) (P±Af, 

where the components of the layer pseudospin are given 
by 



V2 - V4, 



1 A, 



2A(i)' 



(43) 



(P±,l)' = {P.,lf + {Pysf^\{p2A)\\ (44) 

with the convention that pseudospin up is state |2) while 
pseudospin down is state |4) . 



1^1 

{P2a) 



1^3 



1, 



(37) 



-4=2, 

(P4,2) = 2' 



The ground state is an XY layer-pseudospin ferromagnet 
with orbital character n — 1. The Hamiltonian is invari- 
ant with respect to the orientation of the pseudospins 
in the x — y plane so that this phase supports a Gold- 
stone mode. The choice of phase in Eq. ([57]) has the 
pseudospins pointing along the x axis. At finite bias 
Ab < A^\ the pseudospins are pushed out the x — y 



plane i.e. the two layers have unequal population, 
occupation of the four states are in this case 



and 



1^1 ^ 1^3 = 1, 



.(1) 



The 



(38) 



(39) 



B. Inter-orbital-coherent state 

The HFA equations have a separate set of solutions, 
favored at larger bias voltages, in which the lower layer 
is maximally occupied, and the ground state has upper- 
layer orbital coherence instead of inter-layer coherence. 
This solution has 



1^3 = 1^4 = 1, 



U2 = 1 — Ui 



,(2)' 



(45) 



(46) 



and inter-orbital coherence is signalled by the density- 
matrix components 



{P2 



V2 



(47) 



Here 



with inter-layer coherence reflected by 



{P2 



(P4 



V2 - V2- 



We define the critical bias A^^ as the bias at which 
the inter-layer coherence (y02,4) = 0. It is given by 

(1) _ -ix^ (0) + 2£i6 (0)-K2^ 



(40) 



2-4/3 



(41) 



where the Fock interactions Xi (q) and xi (q) are defined 



in Appendix A. As an example, for B 
0.00205 e^Kl 



10 T, A 



(1) 



(2) _ Xi (0) 



2/3 



(48) 



is the critical bias above which all charges in the upper 
layer are transferred to state |2) and the orbital coher- 
ence is lost. In a pseudospin model with the convention: 
pseudospin up for state 1 1) and pseudospin down for state 
1 2), the orbital pseudospin components are given by 



{Sz,k) 

(S^,i> 

{S+,k) 



Vl - V2 

2 ' 
(5j;,i)x+ (S'.y,i)y, 



(pl 



(49) 

(50) 
(51) 
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This orbital-coherent phase has ah pseudospins tilted 
slightly away from the z axis by an angle 



Mixed state 



cos6'i3 = 1-2 



A, 



,(2)- 



(52) 



TT. This critical bias A 



(2) 



At the critical bias A^ , Qb ^ 
is very large; at i? = 10 T, A^-* « he^ j which is near 
the limit of validity of the effective two-band model i.e. 
fiwc < 7i. For Ab > Ag , all electrons are in state |2) 
and there is no further change with bias of the ground 
state. 



by 



E' 



(2) 

HF 



N 



- (1 - /3) ^ + /SAs 
' [xi (0) + xi6 (0) + 2x7 (0)] 



12l 



(53) 



24 

(0)-:ei6 (0)] {S,,k) 



-^[xi {0) + xiG (0)-4x7 (0)] 

-ix4(0) {S^^Kf. 



This energy is independent of the azimuthal angle (p of 
the pseudospin vector. We can thus, without loss of gen- 
erality, take (pi,2) as real. 



40 



Interlayer 
coherence 
in n—1 



30 - 



B ,num 



(A 

« 20 




10 - 



Orbital coherence 
in upper layer 



0.005 



0.015 



A,/(e==/Kl) 



FIG. 3: Phase diagram of the 2DEG in a graphene bilayer 
a.t V = —1. 



Eqs. dUl) and ^ give E 



(2) 
HF 



bias A 



(3) 



HF 
.(3) 



< E'^j^'p above a critical 



.(1) 



a mixed 



The orbital-coherent state (OCS) has an energy given ^(3) 



For biases A e 

state with both inter-orbital and inter-layer coherence 
would be lower in energy than a state with only inter- 
layer coherence. Solving the full Hartree-Fock equations, 
we find that the crossover from the inter-layer coherent 
state to the inter-orbital coherent state occurs contin- 
uously via an intermediate state with both orders. The 
boundaries of the intermediate phase must be determined 
numerically. We find that the boundaries of this mixed 
state (MS) are given on the left by a new critical bias 



< A^-* and on the right by A^-*. 



The intermedi- 



B^num 

ate phase region broadens with magnetic field as can be 
seen in Fig. [S] 

Figure 0] shows the evolution of the inter-layer coher- 
ence (p2,4) and the orbital coherence (pi, 2) with bias at 
magnetic field i? = 10 T. The orbital coherence sets in 
before the inter-layer coherence decreases to zero thus 
creating the mixed state region identified in this figure by 
a non zero value of the density-matrix component (pi,4). 
(Note that (^1,2) is not given by Eq. (|T7|) in the mixed 
state.) The coherences (pi,4) and (^2,3) , which involve a 
mixing of valley and as well as orbital indices, are non- 
zero only in the intermediate mixed-state region of the 
phase diagram. All density-matrix components vary con- 
tinuously with inter-layer bias. 

As can be seen from Figs. [3]and[ll the orbital-coherent 

phase starts at A^' with a finite orbital coherence (^1,2)- 
Were it not for the presence of the inter-layer-coherent 
and mixed states, {pi,2) would start at zero bias and be 
given by Eq. (|47l) for all biaises. The mixed and inter- 
layer coherent states are confined to relatively small inter- 
layer bias voltages; at larger values of As the ground 
state is a relatively simple state with only orbital coher- 
ence. The exploration of collective excitation properties 
of this state is one key objective of this paper. 



VII. COLLECTIVE MODES IN THE 
INTER-LAYER COHERENT STATE 

The density-matrix equations of motions which de- 
scribe the three — —1 spin-diagonal dispersive collec- 
tive modes (see Fig. [2]) are usually simplest when writ- 
ten in the basis of the bonding and antibonding single- 
particle Hartree-Fock eigenstates. For the interlayer co- 
herent state (ICS) (the ground state for A^ < A^^^^^^) 
we find that, 



(I {lu + iS) - Fi) 



( 


ie-"^'ipABi.BQ 








ie'^'^PBo.ABi 









PAB1,B1 









Pbi,abi 




" 




ie'^^PABO^ABl 







\ 


ie-^'^'^PABi.ABo 


) 





(54) 
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where Oq is the angle between the two-dimensional 
wavevector q and the x axis, and B and AB refer to 
the states 



\B,n,X) 
\AB,n,X) 

where n = 0, 1, 



and 



g^\K,n,X)+g+\K',n,X), (55) 
g+\K,n,X)-g^\K',n,X), (56) 



9± 



.(1) 



(57) 



(58) 



FIG. 4: (Color online) Evolution of the inter-layer (p2,4) and 
inter-orbital (pi, 2) coherences (left axis) with bias at filling 
factor V = —1 and magnetic field _B = 10 T. On the right 
axis, one of the non zero coherences (pi, 4) in the mixed state. 



The matrix Fi depends only on the modulus of the 
wavevector q so that the dispersions are isotropic. This 
matrix is given by 



/-A-a^B fH 
-fH A + a^B 
-N + a^J -fj 
fj N-a^J 



-N + a^J -fJ 

fJ N-a^J 
-C - n^D fM 
-fM C + a^D 



-E-a^F K-(T^H 
-K + a^H E + a'^F / 



where 

/ = 1 - a2 (60) 

The variables in this matrix are defined in Appendix B. 

The GRPA dispersions for the three collective modes 
are shown in Fig. [5] for zero bias and a magnetic field 
of B = 10 T. In the limit q ^ 00, the frequencies of the 
dispersive modes correspond to transitions between the 
HFA energy levels indicated in Fig. [2] as expected. Mode 
2 in Fig. [5] is a Goldstone mode consisting of a preces- 
sion of the inter-layer pseudospin Pi around the x axis. 
We refer to it as the inter-layer pseudospin mode (IPM) . 
Mode 1 is an orbital pseudospin mode (0PM) consisting 
of a precession of the orbital pseudospins around their 
local equilibrium position. Mode 3 involves both a layer 
and orbital pseudospin fiip and has a large gap. At fi- 
nite bias, the 3 dispersive modes in Fig. [5] are coupled 
together while at zero bias, modes 2 and 3 completely 
decouple from mode 1 as is clear from Eq. (j59p . 



From Eq. ((5^ . it is easy to see that the dispersion of 
the 0PM at zero bias (given by the 2x2 block in the 
lower right of the matrix Fi) is given by 



o^oPM (q) - 2 (q)-^' (q), (61) 

and has a gap given by (see Appendix B) 

iooPM (0) = i [xi (0) - xi (0) - xi6 (0) + lie (0)] • (62) 

This gap is small but visible in Fig. [S) We find numer- 
ically, that as the bias is increased, the gap loopm (0) 
decreases until it reaches zero at the phase boundary of 
the mixed state. 

The dispersion of the IPM is linear in q for q£ < d/£ 
(with d/i « 0.04 for B = lOT) as in a semiconductor 
bilayer— . This findings differ qualitatively from the = 
—3, case for which we foundii an IPM with q^ dispersion 
and a gapless 0PM. At v = —3 and zero bias, the 15, 0) 
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FIG. 5: (Color online) Dispersion relations of the collective 
modes at = — 1 in the inter-layer coherent phase at zero 
bias and for B = 10 T. The numbers refer to the transitions 
indicated in Fig. [2] 



level is filled. The IPM dispersion in that case occurs 
because the possibility of mixing n = 1 wavefunctions 
with 1^5*, 0) wavefunctions in excited states allows the 
inter- layer phase stiffness to vanish. At = —1, it is not 
possible to make the corresponding admixture. 

In pseudospin language, a finite bias pushes the layer 
pseudospins Pi out of the x—y plane but the Hamiltonian 
of the system remains independent of the orientation of 
the perpendicular component of these pseudospins in the 
x—y plane. The IPM therefore remains gapless for < 

(3) 

num- It acquires a gap in the mixed and OCS. At 
filling factor v = —3, the inter-layer-pseudospin mode 
becomes unstable at finite bias. This indicates that the 
uniform inter-layer-coherent state cannot in fact be the 
ground state at finite bias. We see no such instability at 
ly = -1. 

In the absence of a bias, Eq. (|59p shows that modes 
2 and 3 are coupled through N and J. These interac- 
tk)ns involves the Coulomb interaction matrix elements 
-'^0,1,1,1 i^) -f^o, 1,1,1 (q) (see Appendix A for their def- 
initions). These interactions do not conserve total n = 
or n = 1 quantum numbers. Such interactions do not oc- 
cur in usual semiconductor 2DEG where spin and layer 
pseudospin indices are conserved. 



VIII. COLLECTIVE MODES IN THE ORBITAL 
COHERENT STATE 




FIG. 6: (Color online) Long-wavelength (small q) dispersion 
relation of the collective modes in the inter-layer coherent 
phase at V — —1 for zero bias and for B = 10 T. Mode 1 is 
the orbital pseudospin mode while mode 2 is the inter-layer 
pseudospin mode. 



In this section, we consider collective excitations of 
the orbital coherent state (OCS) which is the ground 
state in the region A^^^ < As < A^^', which covers 
the large region of bias voltages from small values to the 
largest values for which the two-band effective model ap- 
plies. The occurrence of the interesting OCS state at 
high bias voltages is a consequence of competition be- 
tween single-particle and interaction effects as explained 
earlier. By studying its collective interactions we reveal a 
Dzyaloshinskii-Moriya (DM) interaction between orbital 
pseudospins and demonstrate that for large bias voltages 
it drives an instability to an orbital pseudospin spiral 
state. 



A. Electric dipole density 

The fact that (pi,2) 7^ in the OCS implies that there 
is a finite density of electric dipoles in this phase as first 
pointed out in Ref. 12j. To show this, we write the total 
electronic density (including the two valleys) as 

n (q) = UK (q) + riK' (q) , (63) 

where 

riKicO^N^ ^ K,,,{-q)pff {cO, (64) 
j,i=o,i 
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with p^'-^ = /3i,2- The functions Kij (q) are defined in 
Eqs. (I16m9p . In our pseudospin langage for the orbital 
states, we have the relations (for the K valley) 



p"" (q) 


K.K / N , K,K 1 N 

= Po,o (q) + Pi,l (q) 






(65) 


pf (q) 


1 

2 


K,K 1 N K.K 

Po.h (q)-Pi,i 


(q) 




(66) 


p+ (q) 


K.K 1 \ K 1 \ 

^ Po,i (q) = Px (q) ^ 


-^P^ (q) 


(67) 



and so the density operator in the K valley can be written 



as 



UK (q) = N^[l 



(68) 



pk (q) 
) pk,z (q) 

-N^V2iq^epK^^ (q) 
+N^V2iqye -pK y (q) , 

where we have defined (q) — exp (— (7^i'^/4) (q) 
and similarly for ^ with i = x,y, z. 

Now, if the 2DEG is in an external electric field E (r) = 
— V0 (r) , we have for the coupling Hamiltonian 



e / dr n{v) (j) (r) 



(69) 



p,.z (-q) 



iq-r 



1 — I /9, (-q) 



- V2i ((7:,£ Pj- ^ (-q) + qyl p^.^ (-q)) (q) . 

With the electric field in the plane of the 2DEG, this 
coupling can be written, in real space, as 



eN^ ^ / dr[p^.(r)(/.(r) 

j=K,K' 



(70) 



{-p,{v)e~2-p^^^ {v)e) (V.E(r)) 



+y2£ j dr [p^. , (r) (r) - p^,^ (r) Ey (r)] . 

The inter-orbital coherence is zero in the K' valley and 
so {pk',x) j {PK',y) — 0- Eq. (f70| implies that the dipole 
density in the 2D plane is 

(d;, (r)) = -eV2m^ ((p^,. (r)) x - (p^,, (r)) y) 



X ((Pif,^; (q)) X - (pif.y (q)) y) 



(71) 
(72) 



The orientation of the dipole density is set by the phase 



of the density-matrix component (po'-^ (r) ) = (pi^2 (r)) 



which specifies the phase of the spontaneously established 
coherence between n = and n = 1 orbitals in the 
ground state. For our choice of the spontaneously es- 
tablished phase of (pi,2 (r)) in Eq. (H71) . the dipoles are 
oriented along the x axis. Because the K' valley (bot- 
tom layer) Landau levels are maximally filled, there is no 
inter-orbital coherence and therefore no contribution to 
the electric-dipole density from the K' valley. 



B. Effective pseudospin model 

Collective modes dispersions for the OCS are plotted 
in Fig. [71 For relatively small values of they are very 
similar to those represented in Fig. [S] for the ICS; the 
main changes occur at small wavevector as can be seen 
in the inset of Fig. [T] The inter-layer pseudospin mode 
is gapped in the OCS while the orbital pseudospin mode 
is gapless, with a very anisotropic dispersion as shown 
in Fig. [51 We now discuss the physics of the orbital 
pseudospin mode. 




FIG. 7: (Color online) Dispersion relations of the collec- 
tive modes along Qx in the orbital-coherent state at filling 
factor u = —1, bias Ab = 0.0022 / k£ and magnetic field 
B = 10 T. Mode 1 is the Goldstone mode due to orbital coher- 
ence. Modes 2 and 3 involve an inter-layer transition and are 
gapped. The inset shows the dispersion of mode 1 (gapless) 
and mode 2 (gapped) at small QxI. 



At finite wavevector q, the orbital pseudospin mode 
(0PM) corresponds to a precession of the orbital pseu- 
dospins around their equilibrium orientation in the 
ground state as illustrated in Fig. [HI The ground state 
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is described by spinors 



\K, B, X) ^ cos ( -fj \K, 0, X) + sin (^^ j \K, 1, X) , 

(73) 



where 9b G [0,7r] has been defined in Eq. (|52 
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FIG. 8: (Color online) Dispersion relations of the collective 
modes in the orbital-coherent phase at filling factor v = ~1, 
bias As = 0.0022 / and magnetic field B = 10 T. Curve 
(1) qy = 0; curve (2) = qy; and curve (3) q^ — 0. 



The collective mode corresponds to spatially coherent 
rotations of this spinor around its ground state value. 
For this reason, it is convenient to choose the orbital 
pscudospin quantization axis along the direction \B,Q). 
In order to do this, we define 'bonding' and 'antibonding' 
electron creation operators by 



K,K,x = cos 



"ab.k.x 



2 ' ' 



2 

cos I — I c 



,t 

i,x- 



(74) 
(75) 



Below we use the convention that pseudospin up cor- 
responds to state \B,K) and pseudospin down to state 
\AB, K) and we denote the pseudospin by S. 

In our GRPA system of equations for the collective 
modes, the orbital-pseudospin wave mode is decoupled 
from all other modes. Below we follow one possible strat- 
egy for explaining the physics of this mode by comparing 
our microscopic GRPA equation of motion to the equa- 
tions of motion of an effective orbital pseudospin model, 
and using the comparison to identify the effective pseu- 
dospin interactions. Since collective modes correspond to 




FIG. 9: Precession of the orbital pseudospin S around its 
equilibrium position with polar angle 6b and azimuthal angle 
i^B = 0. This orientation is choosen as the orbital pseudospin 
quantization axis in the pseudospin wave model. The angle 6 
is the polar angle measured from the equilibrium orientation 
of the pseudospins. 



small oscillations of the pseudospin around its quantiza- 
tion direction, we can use an effective model which has in- 
teractions only between transverse spins. (Quantization 
direction interactions can be represented as transverse in- 
teractions because of the spin-magnitude constraint.) We 
write the pseudospin effective Hamiltonian in momentum 
(q) space: 



(76) 



where i,j = x^y. Since Jij (r,r') in the real space ver- 
sion of Eq. ([75)1 depends on r — r' only, it follows 
that we can always write Jij{v) = Jji{—v) and hence 
J'ij (q) = Jji (~'3)- Because the real-space interac- 
tions must be real we also have the usual property that 
Jij (q) = Jij (~q)- Combining these two identities we 
can conclude that Jxx (q) and Jyy (q) are real and even 
in q, while Jxy (q) has even real and odd imaginary con- 
tributions. The real parts of Jxy (q) and Jyx (q) are iden- 
tical while their imaginary parts differ in sign. As we 
emphasize further below, the DM interaction is captured 
by the imaginary part of Jxy (q) . 

Using these properties and the commutation relation, 

7V^[5,(q),5,(q')] = i*<5q,-q., (77) 
the equations of motion of the pseudospin model are: 



-iJ. 



xy 



'iJ, 



iJxx 



ijy 



yy 



(78) 



with the dispersion relations 



2ciJ-|- — Z i^ Jxy '^yx) i y ^ Jxx Jyy i^ Jxy ^~ Jyx^ • ('^9) 
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Note that the first term on the right hand side of Eq. 
([75)1 is real and odd in q and that it represents the con- 
tribution of the DM interaction to the collective mode 
frequency. 

Comparing with our microscopic GRPA results for the 
equations of motion and collective mode frequencies, we 
obtain the following expressions for the pseudospin effec- 
tive interactions: 



Jxx 

J xy 
J yx 



+ [A2] - sm\eB) + [A2] )(80) 

-^-svcL(eB)Ai - cos(6is)3 [A2] , 



with 



A2{q) 
A3{q) 

AM) 



hi{q) - xi{q) + Xieiq) ^ hie{q), 
e''"- (he (q) - ^6 (q) 
2hi{q) - 2hi{q) - 2xi{q) + Xi (0) , 
-23* [^e^^■' [q) + hs (q) + X2 (q) - Xs (<?) 



where 9q is the angle between the wavevector q and the 
X axis. All interactions are defined in Appendix A. All 
Hartree and Fock interaction terms in Eqs. (|8ip . hi, hi 
and real and depend only on the modulus of 

q- 

We see from the structure of Eqs. ([50]) that the dis- 



persion relation has the symmetry ll!± ( A^-* — 

7B (As) (from 



uj± (Ab) . Because 6*5 (^A^^^ - Asj = n-l 

Eq. (|52p). it follows that the dispersion at small bias is 
the same as the dispersion near the critical bias as first 
pointed out in Ref. Il2l . 

The physical content of the various terms in Eqs. 
(l80t is most easily identified from their long-wavelength 
forms. From Eqs. ([50)1 wc find that at small q and small 



Jxx (q) 

Jyy (q) 

Jxy (q) 
Jyx (q) 



2/3AB + g^cos2 (6^) 



qism^ (6'q) 



q^£^ (5-6cos2 (6'q)), 




Vf^V V 2^^-^ 9^sin(e'q) 
-qi? sin (6'q) cos (6'q) , 



^^yyl/JAs qi sin (6^) 
—q£ sin (6'q) cos (6'q) . 



(82) 
(83) 

(84) 
(85) 



The pseudospin rotations which change Sx correspond to 
changes in the angle 9 on the orbital pseudospin Bloch 
sphere relative to the ground state value 9b as illustrated 
on Fig. ini For this reason J^x remains finite (is massive) 
as (7 — > when the potential bias is finite unlike the 
other couplings. The terms proportional to qicos^ (0q) 
in Jxx, qisiia^ {9q) in Jyy and g£ sin (6'q) cos (6'q) in J^y 
and Jyx are simply electrostatic interactions between 
changes generated when the dipole orientation varies in 
space. (Recall that the charge density is equal to the 
divergence of the dipole density.) These terms are the 
long-wavelength limits of the Hartree interactions cap- 
tured by the GRPA theory. The imaginary contribution 
to Jxy is the DM interaction whose physics we discuss 
below. The eigenvector for the pseudospin motion, at 
small q and small A^, has Sx/Sy — iy^2q£/ (3Ab sm {9q) 



if 0q 7^ and Sx/Sy = ig^y V27f/(16^As) if 6iq = so 
that the long wavelength collective modes are elliptical 
precessions with minor axis along the massive x direction 
and major axis along the y direction which contributes 
dipolar electrostatic energy. The long wavelength Gold- 
stone collective mode energy therefore has unusual square 
root dispersion: 



(q) = y^2l3ABqi sin(6'q). 
For sin(0q) = 0, we have the linear dispersion: 



a;(q) = ^\/V2^pAi 



(86) 



(87) 



smi9B)~9B^2JAB/A 



We see later that the DM interaction assumes a larger 
importance at larger bias potentials and shorter wave- 
lengths. 

The orbital coherent state occurs at finite bias and 
is preempted at small biases by the interlayer coherent 
state. It is nevertheless interesting to examine the artifi- 
cial limit in which A^ — )■ 0, but layer degrees of freedom 
are still not in play. In that limit all electrons would be 
in the n = orbital, there would be no electric dipoles 
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in the ground state, and the exchange parameters would 
be given by 



Jx 



J, 



3? [^2 



vv 

Jxy 
J yx 



-3[A2]. 



The dispersion relation would then be given by 
^(9)-y^3 (9) -4 1^2 (q)^ 



which is isotropic. The long- wavelength limit dispersion 
would become 



1 /7r\ 1/4 
similar to = — 3 behavior—. 



3/2 



(90) 



to induce an instability of the uniform coherent state. 
When viewed as a classical complex- variable quadratic 
form for Gaussian energy fluctuations, the orbital pseu- 
dospin Hamiltonian, Eq. (|76p. is positive definite pro- 
vided that Jxx{Oi) is positive, Jyyisi) is positive, and 



J xxiOi) J yyiq) > \Jxy{ai 



(92) 



Explicit numerical calculations show that the first two 
stability requirements are always satisfied, but that be- 
cause of the DM interaction, the third is not satisfied 
when sin(0B) is large. In the GRPA we find that the 



(89) 0PM first becomes soft at qyl ~ 2 when As = 



{DM) 



A^^VlO- (We have Ob « 41° at this value of As-) Fig. 
[TOl shows the instability of the orbital-pseudospin mode 
at A^*^-*. We remark that the instability occurs at a 
positive (negative) value of qy in cj+ (w_) . The higher- 
energy collective modes (not shown in the figure) show 
no sign of instability. 



C. Moriya interaction and spiral state instability 

As explained previously the DM interaction is captured 
by the imaginary part of Jxy (q) . When this contribution 
to the orbital pseudospin Hamiltonian is isolated it yields 
an interaction of the standard^'^ DM form: 



Hdm = iN^ 



q 



3 [J: 



xy 



(q)](S(q)xS(-q)).2. (91) 



An examination of Eqs. (|8T|) shows that this interac- 
tion is not due to electrostatic dipole interactions, and 
instead to the exchange vertex corrections i.e. to the 
interactions X2 {q) and xg (g). In Appendix C, we ana- 
lyze the exchange energy of quantum Hall ferromagnets 
quite generally and show that DM interactions are the 
rule rather than the exception when the two-states from 
which the pseudospin is constructed have the same spin. 
The physics of exchange interaction contributions to col- 
lective mode energies is most simply described by mak- 
ing a particle-hole transformation for occupied states, as 
discussed in Appendix C. The exchange interaction at 
momentum p can then be related^^ to the attractive in- 
teraction between an electron and a hole separated by 
. DM interactions occur when the pseudospin state of 
the electron or hole give rise to cyclotron orbit charge dis- 
tributions which do not have inversion symmetry, a prop- 
erty that holds here because of dipole formation. These 
distortions of the cyclotron orbit are irrelevant when p 
is very large but become important for p ^ . In the 
case of a simple parabolic band = I quantum Hall fer- 
romagnet, for example, this picture^^ provides a simple 
understanding of the full spin- wave dispersion. 

(2) 

The DM interaction is strongest at Ab/Aq = 1/2, 
i.e. when s'm{9B) = 1. We have found that over a broad 
range of A^ values the DM interaction is strong enough 




FIG. 10: Dispersion relation in the orbital-coherent state at 
the Dzyaloshinskii-Moriya instability. Filling factor u = —1 
and bias As = 0.51 j kI. 



The eigenvector with positive frequency of Eq. ([78 
has 



Sx 

Sq, 



J X 



J y 



^JxxJyy {.Jxy ^yx) 



(93) 



It follows, using Eq. (|92|) . that, at the DM instability, 
the energy is lowered by forming coupled density-waves 
in X and y pseudospin components with: 



Sxiq) _ Jxyiq) _ Jyyiq) 



Syin) 



Jxx{q) J* (q) 



(94) 
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The real part of the Jxy couphng, due mainly to the 
dipole electrostatic energy, is very small at the instabil- 
ity wavevector because the = Landau level cannot 
support rapid spatial variation, as we have verified by ex- 
plicit calculation. Because Jj,j,(q) is real, it follows that 
Sx and Sy spatial variations are out of phase by nearly 
exactly 7r/2. If the magnitudes of the Sx and Sy com- 
ponents were identical, this would imply a spiral ground 
state. Because Jxx(q) and Jyy{ci) are not identical at 
the instability, the spiral is somewhat distorted. It must 
be kept in mind, however, that the DM instability may 
be preempted by a first order transition to a state with 
lower energy and a more complex pseudospin pattern. A 
fuller exploration of the properties of these states, includ- 
ing their properties in the presence of an external electric 
field, is beyond the scope of the present work. 



IX. COLLECTIVE MODES IN THE MIXED 
STATE 

The mixed state occurs between the inter-layer coher- 
ent and inter-orbital coherent phases as shown in Fig. |3l 
The width of this region in the phase diagram increases 
with magnetic field. In this phase, all order parameters 
(pi.j) are finite so that this phase has both inter-layer and 
inter-orbital coherence. We show the behavior of some 
of these order parameters with bias in Fig. |4l The order 
parameters (^1,4) and (^2,3) which flip both valley and 
orbital indices are non zero in this phase. 

The collective modes in the MS are obtained numeri- 
cally by solving Eq. (p6| . The MS has three dispersive 
modes, as in the other two phases we studied. The dis- 
persions of these modes differ from the dispersions in the 
other two phases at small wavevector only. We show in 
Fig. fTlT a) the dispersion of the inter-orbital mode and 
in Fig. Illf b) the dispersion of the inter-layer coherent 
mode along the x or y axis and at 45° {xy) from the x 
axis. In contrast with the other two phases, both modes 
are now gapped for not at the boundaries of this 
phase. As in the inter-orbital phase, the dispersions are 
highly anisotropic. We show in the next section that this 
phase has a distinct signature in the microwave absorp- 
tion spectrum. 



X. MICROWAVE ABSORPTION 

The collective modes discussed in the previous sections 
can be detected in microwave absorption experiments, as 
we now show. We write the current operator, projected 
onto = and valley K, as 



dAf 



(95) 
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FIG. 11: Dispersion relation of the (a) inter-orbital and 
(b) interlayer coherent modes in the mixed phase at As = 
0.00195 c^/kI and _B = 10 T along directions x,y and xy 
(45° from the x axis). 



second quantization, the total current is given by 

J^K = J drM'y (r) j^A'^-fK (r) , (96) 

with the field operators defined in Eqs. (|7l8p and ^ = ±1. 
Wc find that 

r- el 

3^K - e^/3AB-A^^ {p^K.y (0) X + p^K^x (0) y) , (97) 
where 



(98) 
(99) 



PK,.(0) = \ [p^f{Q) + p'^f{Q) 
PK,yiO) = ^(poT(0)-P^!f (0)) 



and similarly for px' ,x/y (0) ■ The same result for J^k can 
be obtained by calculating the polarization current 



where Af is the vector potential of the external electro- 
magnetic field, i = x,y and is given in Eq. In 



J«if = ^djK (q) , 



(100) 
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with the dipole density d^K (q) defined in Eq. (|72)) and 
using the Heisenberg equation of motion —iHd/dt {...) — 



(...) 



We note that the orbitally coherent state 

has spontaneous currents in its ground state, a very ex- 
ceptional property. It appears hkcly that, in finite sys- 
tems, these currents should flow perpendicular to sys- 
tem boundaries, forcing domain structures in the pseu- 
dospin magnetization texture, consistent with expecta- 
tions based on the electrostatic energy of the dipole den- 
sity. 

We define the total current-current correlation func- 
tion Matsubara Green's functions as 



-|(TJ„(r)J^(0)) 



(101) 



71 / 47rm*2^4 



i,j=K,K' 



where Ja — Jk.ol + JK',a and 



{t)^-{Tp,,^ (0,r)p,-^ (0,0)) 



(102) 



with i,j = K,K' and a,f3 = x,y. Note that a, /? are 
defined so that x — y,y = x. 

The microwave absorption for an electric field oriented 
along the direction a is given by 



Pa(w) = 



ic 

h 
1 

2 



xl.j„ (^) 

71 



(103) 



E 

i,j=K,K' 



R 



i5 



E, 



0' 



where we have assumed a uniform electric field E (t) — 
i^oSe""^* and taken the analytic continuation ifln — > 
u! + i5 oi XJc,Jc (0, i^n) to get the retarded response func- 
tion. The response functions Xpi^p^^ are calculated in 
units of h/ {e^ / ntj so that (w) is the power absorbed 
per unit area. In Eq. (jl03p we have neglected a diamag- 
netic contribution to the current response which becomes 
important at low frequencies. 

Our GRPA correlation functions are given by Eq. 
and numerical results for the absorption in the inter-layer 
coherent phase are shown in Fig. [12] Exactly the same 
result is obtained, in this phase, if the electric field is set 
in the y direction, i.e. the absorption is isotropic. We 
see that the signal in the absorption is at a frequency 
corresponding to the orbital-pseudospin mode (see Fig. 
|6]). The frequency v of this mode at q = decreases 
with bias while the absorption intensity increases with 
bias. Since e'^/hnl - 2.7 x 10 (T) GHz (using k = 5 
for graphene on Si02 substrate), the frequency of the 
orbital pseudospin mode at zero bias \s v k, 3.4 GHz i.e. 
in the microwave regime. Note that mode 3 in Fig. [S] 
is also present in the absorption at finite bias and that 



its frequency has a higher value outside of the microwave 
regime. 

The absorption in the mixed state is shown in Fig. 
[T^ on a logarithmic scale. The three lower peaks show 
the absorption from the inter-layer coherent mode which 
is gapped in the mixed state. The other three peaks are 
from the inter-orbital mode. The gaps in these two modes 
increase with bias until they reach a maximum around 
As = 0.00195 for B 10 T. The gaps then decreases 
with bias. The intensity of the absorption increases with 
bias for both modes. In Fig. [131 the electric field is 
set along the y axis. The absorption is at least 100 times 
lower if the electric field is oriented along the x axis i.e. it 
is highly anisotropic in the mixed state. With our choice 
of phase for the ground state of the mixed state, the 
pseudospins (and so the electric dipoles according to Eq. 
(1711) ) are oriented in the y direction. Their motion for 
q = is an oscillation in the x — y plane about the y axis. 
For such a configuration, the electromagnetic absorption 
is strongest for fields along the y axis, which is what we 
observe in our calculation. 

In the orbital-coherent phase, the orbital-pseudospin 
mode is gapless and decoupled from the two other gapped 
modes. Since the orbital-coherent mode couples strongly 
to external electric fields, we can expect anomalous low- 
frequency absorption in this state, similar to the Drude 
absorption of a metal. This interesting and unusual ab- 
sorption feature is likely to be highly sensitive to dis- 
order. Its detailed analysis lies beyond the scope of 

(2') 

the present paper. Above , the ground state has 
V2 = = = 1 and there is no orbital coherence any- 
more. The 0PM then has a gap that is proportional to 

Ab — A^ . The 0PM becomes visible in the absorption 
in this phase while the other modes do not. 

In summary, we see that each of the three phases in 
the phase diagram at v = —1 has a different signature in 
the microwave absorption spectrum. We remark that the 
frequencies in Fig. [T3|are quite small. But, they can be 
increased by increasing the magnetic field. For example. 
Fig. [T4l shows the gap in the inter-layer and inter-orbital 
pseudospin modes in the three phases at a magnetic field 
of 40 T. 



XI. CONCLUSION 

We have studied the phase diagram and collective ex- 
citations of a spin-polarized bilayer graphene 2DEG at 
filling factors v = —1 and = 3, as a function of a bias 
electric potential which shifts electrons between layers. 
Our study is based on the Hartree-Fock approximation 
for the mean-field ground state, and the GRPA for the 
collective modes and response functions. 

We predict phase transitions between the following se- 
quence of states with increasing potential bias: (1) an 
inter-layer-coherent state (ICS) with a zero gap inter- 
layer pseudospin mode (IPM) and an orbital pseudospin 
mode (0PM) with a small gap. This gap can be detected 
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FIG. 12: (Color online) Microwave absorption from the or- 
bital pseudospin mode in the inter-layer phase at different 
values of the bias. The absorption is zero at zero bias. The 
second(third) peak has been multiplied by 10(100). 



FIG. 13: (Color online) Absorption spectrum in the mixed 
state at _B = 10 T for an electric field in the y direction and 
for biases As = 0.00185, 0.00195, 0.0020 in units of e'^/ni 



in microwave absorption experiment. Its frequency de- 
creases with bias while its intensity increases with bias; 
(2) a mixed state (MS) with both inter-layer and inter- 
orbital coherence. Both the IPM and 0PM are gapped 
and visible in microwaves in this phase. Moreover, the 
intensity of the absorption is highly sensitive to the direc- 
tion of the external electric field; (3) an orbital coherent 
state (OCS) with orbital coherence concentrated in one 
layer only. (The second layer is completely filled.) The 
OCS state is a very simple one in which the low poten- 
tial layer is filled and the high potential layer has a gap 
induced between it's two Landau levels by spontaneously 
establishing coherence between states with n = 1 and 
n = orbital character. This state has a number of 
quite unusual properties, including electric dipoles and 
associated uniform currents in its ground state. The 
OCS has a gapless (Goldstone) 0PM and a gapped IPM. 
Both modes are absent from the absorption spectrum. 
The dispersions of the collective modes are also highly 
anisotropic in this phase. The phase is unstable at a fi- 
nite wavevector due to the presence of a Dzyaloshinksii- 
Moriya exchange interaction. We believe that the insta- 
bility will lead to the formation of a ground state with a 
non-uniform pseudosopin pattern. 

These properties are associated with competition be- 
tween an electron-electron interaction term which favors 
n = orbital occupation, and single-particle terms in the 
effective two-band model of Ref. fsl which favors n — 1 
orbital occupation more strongly at larger inter-layer po- 
tential difference. We note that this band model ne- 
glect the 74 hopping term, which connects sites Ai(Bi) 
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FIG. 14: (Color online) Gaps in the inter-layer pseudospin 
mode (IPM) and in the orbital pseudospin mode (0PM) as a 
function of bias ai v = —1 and B = 40 T. 



and A2{B2)^- Strictly speaking, the results obtained 
in our paper are valid if 74 — 0. The effect of the 74 
term is to add a negative correction to the energy of the 
n = 1 orbital which is independent of the valley index, 
magnetic field and bias. More precisely: Ej^ n=i x = 



/3Ab + C,EK'.n=l,X 



\Ab + 13Ab + C with 
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C = -\/6^74ao/^. It is not clear from the literature what 
is the precise value of the 74 parameter. In Ref. [l^, 
the value 74 = —0.12 is obtained by comparing the tight- 
binding model for bilayer graphene with the Slonczewski- 
Weiss-McClure tight-binding model for bulk graphite^*^. 
In bulk graphite, 74 w 0.044 eV. With 74 = -0.12, the 
correction C < for Ab > 0.02 e'^/ni. It follows that 
this term may have an effect on the existence of the inter- 
layer coherent and mixed phases, and certainly on the lo- 
cation of their phase boundaries, A^^^^^ and Ag \ since 
level n = 1 may be below level n = in both layers. The 
Dzyaloshinksii-Moriya transition, however, occurs at a 
large bias of A^ « 0.5 e^/ni ma region where the single- 
particle energy of level n = 1 is already well below that of 
level n = 0. For this reason we believe that the physics in 



this region should not be affected by 74. More numerical 
work is needed, however, to assess the precise influence 
of this term. 
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Appendix A: HARTREE AND FOCK 
INTERACTIONS 



We first give the definitions of the Hartree, H, and 
Fock, X, interactions in Eqs. (I12I13I14I15P : 



hi{q) = iJo,o,o,o (g) = — ;A (g) 



(Al) 



h2 (q) = i/0,0,0,1 (q) = -^^''"^ (9) ' (A2) 



/14 (q) 


= ^0,04,1 (9) = 


1 / q^P\ 


AW, (A3) 


K (q) 


= -^0,1,0,1 (q) = 


\qie^^'-K{q), 


(A4) 




^0,1,1,0 (9) = 


^'ZM(g), 


(A5) 


hs. (q) 


= -ffo,i4,i (q) = 


)A(g), 


(A6) 



h^, (q) = i?i,M,i (9) = ^ (1 - ^ ) A (9) (A7) 



and 



h3 - i/oAi,o(q) =-/i2, (A8) 

hs = i/o,i,o,o(q) = -h2, (A9) 

h9 = H,^o,o,o{(i)=K, (AlO) 

hio = Hi^o.o,! (g) = (All) 

hii = ifi,o,i,o(q) =/i5, (A12) 



19 



hi2 = Hi,o,i.i (q) = -h;, (A13) 
hi3 = -ffi, 1,0,0 (q) = hi, (A14) 
hi4 = 1,0,1 (q) = ~hs, (A15) 
/»15 = i?i,i,i,o(q) = /i8: (A16) 
where 9^ is the angle between the wavevector q and 

— ^ — 1 . The interactions 

hn (q) are obtained by multiplying ft,„ by e^'^'' where 
d is the inter-layer separation. The interactions hn (q) 

and hn are obtained by removing the term i and the 
phase factor e^'^^ or e^^*^"! in hn and For example 
h2 - -A (g) /V2 while /i2 (q) = -^e^'^'^A (g) . 
The Fock interactions are defined by 

XI (q) = Ao,o.o,o (q) = / dye-^'/Vo ((?£y) , (A17) 
Jo 



X2 (q) = Xo,o,o,i (q) = ^e^'" dyye-y"/\h {qty) , 

(A18) 

/■oo / 2 \ 

X4 {q) = ^0,0,1,1 (q) = y ( 1 ~ y ) e"^'/'-/o (g^y) , 

(A19) 

xe (q) = Xo,i,o,i (q) = 2^'''" / dyy^e-^'^Va , 

° (A20) 

xr (q) = Xo,i,i,o (q) = 2 / dyy\-y"^^Jo {qly) , 



(A21) 



a;8(q) = Xo,ia,i (q) = -^e'''- dyy (A22) 



poo / 2 \ ^ 

xi6 (<?) = Xi,ia,i (q) = J dy(l-^\ e-«'/Vo (g^) , 
° (A23) 

and 



a;3 


— -'^o, 0,1,0 (q) 




(A24) 


X5 


— -'^o, 1,0,0 (q) 


= -X2, 


(A25) 


Xg 


— -'^i, 0,0,0 (q) 




(A26) 


XlO 


= -^1,0,0,1 (q) 


= XT, 


(A27) 


Xii 


= ^1,0,1,0 (q) 




(A28) 


Xl2 


= ^1,0,1,1 (q) 


— Xg, 


(A29) 


Xl3 


= -'^i, 1,0,0 (g) 


= X4, 


(A30) 


Xi4 


= -'i^i, 1,0,1 (q) 


-Xs, 


(A31) 


Xl5 


= -''^1,1,1,0 (q) 


— — Xg. 


(A32) 



The interactions x^i are obtained by multiplying the inte- 
grand by e~y'^/^. The interactions Xn and are obtained 
by removing the imaginary term i and the phase factor. 
The combinations: 



Hi hi hi , 

Ti hi -\- hi , 

^i — Xi H~ Xi , 

Ui Xi Xi . 



(A33) 
(A34) 
(A35) 
(A36) 



To define Hn,Tn, Xn,Un, we follow the same procedure 

as for hn , hn , Xn T Xn- 

Some useful constants: 



(A37) 



^1 (0) = \/2' 



X4 (0) 



^7 {0)^T;\h 



1 hi 
2V 2' 

1 

2V2' 



2:^16 (0) 



3 hF 
iV 2- 



(A38) 



(A39) 



(A40) 



Appendix B: MATRIX Fi FOR THE 
INTER-LAYER-COHERENT MODES 

The collective modes a.iv = ~l and < A^' involve 
the matrices Fi (q) in Eq. (|54p. The elements of this 
matrix are defined by 

Aiq) = ^4 (0) - 5?i (0) - (0) - (0) (Bl) 
+X4 (q) - Hi (q) + H4 (q) , 

B (q) = 2A| (/3 - 1) + Hi (q) - H4 (q) - U4 (q) (B2) 
+xi (0) - ^xi (0) - X4 (0) + xiG (0) + 2^, 

C (q) = -i/i6 (g) + A16 (g) - 2^16 (0) , (B3) 

D{q) = 2A%{2/3-l) + Hieiq)-UiUq) (B4) 
+25;i6 (0) - -XI (0) + 2-, 

(g) = T4 (q) + X4 (q) - Ti (g) -|- 5?i (0) (B5) 
-X4 (0) - ^xi (0) - xie (0) , 
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Fiq) = 2/?A^ - Fi (g) + F4 (9) + C/4 (g) (B6) 
+X4{0)-xi {Q)-^xi {Q)+xie (0) , 



G{q) = A's-H,{q) + H,{q) + U4{q) (B7) 
+^4(0)-^! (0) + ixi (0)-^, 



Hiq)=He{q)~Ue{q), 

Jiq)^Hsiq)-Us{q), 

K{q)=Xe{q)-fe{q), 

M iq) = H,e (q) - C/16 («) , 

N{q)^Hs{q)~Xs{q). 
To get the gap given by Eq. (15^ . we use 



(B8) 



(B9) 



(BIO) 



(Bll) 



(B12) 



|£;(o)| = 



-Ja/I + (0)-xi6 (0) 



- X, (0) - (0) - x^e (0) + 5?i6 (0) ,(B13) 
K{0) = 0. (B14) 



Appendix C: EXCHANGE ENERGY 
PSEUDOSPIN DEPENDENCE 



It is possible to derive a rather general and instruc- 
tive expression for the exchange energy of a pseudospin- 
1/2 quantum Hall ferromagnet, for the case in which the 
pseudospin texture varies in one direction only. In the 
following we take this direction to be the x direction and 
choose a Landau gauge in which the guiding centers or- 
bits are localized as a function of this coordinate. If we 
are interested only in the dependence of energy on pseu- 
sospin texture we can assume that every guiding center 
orbital is occupied by one electron, but leave the pseu- 
dospin of that orbital arbitrary. It is not necessary to 
immediately specify the orbital character of the states 
that form the pseudospin and we refer to them for the 
moment as state A and state B. In the case of immediate 
interest in this paper, state A has n = orbital character 
and state B has n — 1 orbital character. We discuss some 
other examples below. 

The pseudospin texture can be specified by the two- 
component pseudospinors at each guiding center. 



(CI) 



or by the guiding center dependent direction cosines 
of the pseudospin orientation: {n^,n^^,n^) = 



{s'm{9x) cos{(j>x), s'm{0x) sm{(j>x), cos{9x)) where 
0x and (px are the pseudospin orientation polar and 
azimuthal angles. When a specific gauge choice is conve- 
nient we use — cos(6'/2) and — sin(6'/2) exp{i(f>). 

The exchange energy of a pseudospin quantum Hall 
ferromagnet is 



-Ex 



-- J2 {X,X'\V,,\X',X) 

X,X' 



(C2) 



The dependence of exchange energy on pseudospin tex- 
ture can be exhibited explicitly by using the property 
that each exp(ik • (ri — r2)) term in the Fourier expan- 
sion of the two-particle interaction matrix elements can 
be separated into factors that depend on ri and Y2 inde- 
pendently. It follows that 



X,X' 



>'X 



' '^X'X ■ 



(C3) 



where the Greek indices a = c,x,y, z, n'^ = I represents 
the filling factor of the guiding center states. 



-?y) 



-^xtx = ^E°>^P(-9'^V2) r-'i<D Vi<D 5x'-x, 

(C4) 

i is the magnetic length, V (q) is the Fourier-transform of 
the electron-electron interaction and /"'^(q) is an inter- 
action form factor. Eq. (jC4p follows from the following 
expression for the plane-wave matrix elements: 

(X'|exp(^q•r)|X) - exp(-(7'^V4) 



^ zi:,F^'^(q)zi, (C5) 



I' J 



where the overbar accent denotes complex conjugation. 
The character of the cyclotron orbitals in the two nearly 
degenerate Landau levels is captured by the single- 
particle form factors ''^(q)- For example if pseu- 
dospin / has orbital Landau level index n/, ''^(q) = 
Fnji,ni{f\} where -F'ri'.n(q) is the familiar two-dimensional 
electron gas Landau-level form factor—, commonly used 
in the analysis of many-different properties in the quan- 
tum Hall regime. Note that F^''^{q) = F^'^'(-q) For 
the example of quantum Hall ferromagnetism discussed 
in this paper, ua = and ub — 1 when the bias poten- 
tial is strong enough to yield complete layer polarization. 
(The analysis in this section does not apply when both 
orbital and layer degrees of freedom are in play.) When 
the A and B orbitals have opposite spins, a common oc- 
currence in quantum Hall ferromagnetism, F^'^[c() = 0. 
Examples in which the A and B orbitals are centered in 
different two-dimensional layers require a slight general- 
ization of the present discussion which we do not explic- 
itly address. 

Using Eq. (jC5l) we find (leaving the wavevector depen- 
dence of the single-particle and interaction form factors 
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\FB.B,2 ^ lpA.B,2 ^ lpB.A\2 



implicit) that 

4/0'^ = -2Im[F'^'^^^^^^ + F^^^F^^^] 
4ja,a = 2Re[F^'^F^'^ - F^'^F^'-^] 

^fz.z ^ |^A,A|2 |^S,B|2 _ |^A,B|2 _ |^S,A|2 

4^-2' = 2Ini[F'^'^F^^-^ +F^''^F"^'-^] 
4j^,2 = 2Re[F^'^F-^''^ - F-^'-^F^'-^] 
4ja,2 = _2Im[F-^'^F^'^-F^^-«F^^-«]. (C6) 

These resuhs are obtained after replacing the pseu- 
dospinors by pseudospin direction cosines using 

2z^z^ ^ 1 + n' 
2z^z^ = 

2z^z^ = rf^iny (C7) 



Note that because the exchange energy is real, all the 
interaction form factors are real. The diagonal interac- 
tion form factors are even functions of q, while the off- 
diagonal form factors have both even and odd contribu- 
tions. The interaction form factors capture the influence 
of the shape of the pseudospin state dependent cyclotron 
orbits on exchange energies. In the limit q — 0, or- 
thogonality implies that F^'^ ^ ^ab- If follows that 
/">/3(q = 0) = 5<,,^/2. 

The exchange energy expression can be written in an 
alternate form by Fourier transforming along the direc- 
tion perpendicular to the guiding center orbitals, defining 



1 



E 

X 



exp(— ipX) n\ 



(C8) 



where A^^ = ^-kP^ j A is the number of guiding center or- 
bitals in a Landau level. The exchange energy then takes 
a form from which spin-wave dispersions can be simply 
read off: 



1 



(C9) 



with the momentum space exchange integral given by 



•Jp 



(2^ 



exp(-w^^2) g^p(_^2^2/2)/",^*(q) Via). 

(CIO) 



Because /"'^ is real, J'^p = Jp'^ ■ Similarly, for the 
diagonal elements of J, the property that the diagonal 
elements of / are even functions of q implies that ''^ = 

J-p- Combining both properties, we conclude that the 
diagonal momentum space interactions are real. The off- 
diagonal elements, however, can have both even real and 
odd imaginary contributions. For 3> 1, the integra- 
tion over q in Eq. (jC10[) has contributions from small q 



only, allowing us to set exp(-g^^2/2) /"-^(q) 5afi/2. 
The integral can then be recognized as an inverse Fourier 
transform of the electron-electron interaction from mo- 
mentum space back to real space. It follows that. 



a.fi 



5a 



fit 



2e\p\P' 



(Cll) 



This is the familiar peculiarity of quantum Hall sys- 
tems in which exchange interactions at large momenta 
are most simply understoocti^ after a particle hole trans- 
formation which converts them into Hartree interactions 
between electrons and holes. In strong fields the particle 
and hole in a magnetoexciton with momentum p are sepa- 
rated in real space by pl"^ . For very large p, the separation 
between particle and hole is larger than the cyclotron or- 
bit sizes and the interaction is approximately given by the 
interaction between point charges. For smaller values of 
p the size and shape of the cyclotron orbit, represented in 
momentum space by exp(— g^£^/2) /"'^(q) becomes im- 
portant. 

Imaginary off-diagonal contributions to J!^'^ are re- 
sponsible for DM-like interactions in quantum Hall fer- 
romagnets. In the example discussed in this paper for 
instance, the A orbital has orbital Landau level index 
n = 0, while the B orbital has the same spin and 
n = l. It follows that F^'^ = LQ{q^l^l2) and F^^^ = 
L\{(^P' jl) , where i„ is a Laguerre polynomial, and 
that F^'^ — ig£exp(i0q)/-y/2 where f?q is the orienta- 
tion angle of q. It follows that Z^'** vanishes. The gen- 
eralized random-phase-approximation that we use in the 
main text for collective mode calculations is equivalent 
to linearized pseudospin-wave theory. For fluctuations 
around a ground state with pseudospin orientation polar 
angle Qb (as discussed in the main text), the J^'^ and 
J^^^ pseudospin interactions give rise to DM-like interac- 
tions whose strength is proportional to sin(6's). The DM 
interactions can be viewed in particle-hole language as a 
consequence of the property that the interaction between 
an electron with pseudospin in the x — y plane at X and a 
hole with pseudospin in the z direction at X' is not invari- 
ant under the interchange of X and X' . This property 
reflects the pseudospin dependence of dipole and other 
higher order multipoles in the cyclotron orbit cloud of an 
electron with a particular pseudospin. 

The analysis presented in this section, can be applied 
to any quantum Hall ferromagnet. One elementary ex- 
ample is the case where the pseudospin orbitals share 
Landau level index n and have opposite spins. In this 
case 



/ 



^ Ea,^{L^{qH^ 12)) 



(C12) 



The pseudospin model has isotropic Heisenberg interac- 
tions and the pseudospin-wave excitation energy expres- 
sion ujp = 2(Jp^o ~ Jp)^ follows immediately from an 
expansion of the exchange energy function to quadratic 
order. 



